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. We predict the formation and superfluidity of polaritons in an optical microcavity formed by 

excitons in gapped graphene embedded there and microcavity photons. The Rabi splitting related 
to the creation of an exciton in a graphene layer in the presence of the band gap is obtained. 
The analysis of collective excitations as well as the sound velocity is presented. We show that the 
superfluid density n s and temperature of the Kosterlitz-Thouless phase transition T c are decreasing 
functions of the energy gap. 
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C3 ■ I. INTRODUCTION 

Up today all theoretical and experimental studies have been devoted to Bose coherent effects of two-dimensional 
(2D) excitonic polaritons in a quantum well embedded in a semiconductor microcavity [H-0]. To obtain polaritons, two 
mirrors placed opposite each other in order to form a microcavity, and quantum wells are embedded within the cavity 
| at the antinodes of the confined optical mode. The resonant exciton-photon interaction results in the Rabi splitting 
of the excitation spectrum. Two polariton branches appear in the spectrum due to the resonant exciton-photon 
coupling. The lower polariton branch of the spectrum has a minimum at zero momentum. The effective mass of the 
lower polariton is extremely small. These lower polaritons form a 2D weakly interacting Bose gas. The extremely 
light mass of these bosonic quasiparticles at experimentally achievable excitonic densities, results in a relatively high 
critical temperature for superfluidity, because the 2D thermal de Broglie wavelength, which becomes comparable to 
^ ■ the distance between the bosons, is inversely proportional to the mass of the quasiparticle. 

Recently there were many experimental and theoretical studies devoted to graphene known by unusual properties 
in its band structure [H, Q . Due to the absence of a gap between the conduction and valence bands in graphene, 
the screening effects result in the absence of electronic excitations in graphene. Today we achieved different ways to 
t^J- . obtain a gap in graphene. For example, the gap in graphene structures can be formed due to magnetic field, doping, 
electric field in biased graphene and finite size quantization in graphene nanoribbons. A gap in the electron spectrum 
r**"" , in graphene can be opened by applying the magnetic field, which results in the formation of magnetoexcitons 0]. 

Excitons in graphene can be also formed due to gap opening in the electron and hole spectra in the graphene 
I/"") \ layer by doping [8(. There were a number of papers devoted to the excitonic effects in different graphene-based 
structures. Significant excitonic effects related to strong-electron-hole correlations were observed in graphene by 
measuring its optical conductivity in a broad spectral range [9j. The observed excitonic resonance was explained 
within a phenomenological model as a Fano interference of a strongly coupled excitonic state and a band continuum Q . 
^ ' The electron-hole pair condensation in the two graphene layers have been studied in Refs. fToj - Tl^ . The possibility of 
formation of edge-state excitons in graphene nanoribbons is caused by the appearance of the gap in the electron energy 
spectrum due to the finite size quantization. A first-principles calculation of the optical properties of armchair-edged 
graphene nanoribbons with many-electron and excitonic effects included was presented in Ref. [H. In Ref. [H were 
studied the optical properties of zigzag-edged graphene nanoribbons with the spin interaction. It was found that 
optical response was dominated by magnetic edge-state-derived excitons with large binding energy. First-principles 
calculations of many-electron effects on the optical response of graphene, bilayer graphene, and graphite were described 
in Ref. 1 171 It was found that resonant excitons were formed in these two-dimensional semimetals. The other mechanism 
of electronic excitations in graphene can be achieved in biased graphene, where the energy band gap is forme d by 
applied electric field. A continuously tunable bandgap of up to 250 meV was generated in biased bilayer graphene [18j . 
It was shown that the optical response of this system is dominated by bound excitons [r| . 

According to Ref. a tunable gap in graphene can be induced and controlled by hydrogenation. The excitons 
in gapped graphene can be created by laser pumping. The superfluidity of quasi-two-dimensional dipole excitons in 
double-layer graphene in the presence of band gaps was proposed recently in Ref. [2ll . 

The Bose-Einstein condensation (BEC) of magnetoexcitonic polaritons formed by magnetoexcitons in graphene 
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embedded in a semiconductor microcavity in a high magnetic field in a planar harmonic potential trap was studied 
in Ref. [22T However, the interaction between two direct 2D magnetoexcitons in graphene is negligibly small in a 
strong magnetic field, in analogy to 2D magnetoexcitons in a quantum well [23J]. Therefore, the superfluidity of 
magnetoexcitonic polaritons in graphene in this case is absent, since the superfluidity is caused the sound spectrum 
of Bosc collective excitations due to the exciton-exciton interaction which is negligible in graphene in a high magnetic 
field. 

In this paper we consider the direct 2D excitons formed in a single graphene layer in the presence of the band 
gap and predict the superfluidity of polaritons formed by these excitons and microcavity photons, when the graphene 
layer is embedded into an optical microcavity We obtained the corresponding superfluid density and temperature for 
the Kosterlitz-Thouless phase transition due to the superfluidity of microcavity polaritons. 

The paper is organized in the following way. In Sec. [TT] we presented the Hamiltonian of excitons in a graphene 
layer embedded in an optical microcavity. In Sec. IIIII we obtain the excitonic Hamiltonian which is the sum of the 
Hamiltonian of non-interacting excitons in the gapped graphene and Hamiltonian that describes the exciton-exciton 
interaction. The Hamiltonian of photons in a semiconductor microcavity is given in Sec. IIVI In Sec. |V]the Hamiltonian 
of the harmonic exciton-photon coupling in the gapped graphene is derived and corresponding Rabi splitting constant 
is obtained. The study the condensation of a gas of microcavity polaritons, the density of the superfluid component, 
as well as the Kosterlitz-Thouless temperature are presented in Sec. IVII Finally, the discussion of the results and the 
conclusions follow in Sec. I VIII 

II. HAMILTONIAN OF GAPPED GRAPHENE EXCITONS IN MICROCAVITY 

The total Hamiltonian H to t of the system of 2D excitons in gapped graphene embedded in an optical microcavity 
and 2D microcavity photons can be written as 

Htot = H ex + H p h + H ex -ph , (1) 

where H ex is the Hamiltonian of excitons in graphene in the presence of the gap, H p h is a Hamiltonian of photons in 
a microcavity, and H ex - P h is a Hamiltonian of exciton-photon interaction. 

The Hamiltonian of 2D excitons in the graphene in the presence of a gap is given by 

Hex — H^x Hex — ex: (2) 

where Hex is the Hamiltonian of non-interacting 2D excitons in gapped graphene and H ex -ex is the Hamiltonian of 
the exciton-exciton interaction. 

The Hamiltonian of non-interacting excitons in gapped graphene is given by 

p 

where bp and Sp are excitonic creation and annihilation operators obeying to Bose commutation relations and £ ex (P) 
is the energy dispersion of a single exciton in a graphene layer. 

The Hamiltonian of the exciton-exciton interaction H ex -ex in graphene in the presence of a gap is given by 

Hex-ex = ^ ^P+q^P'-qMP'' (4) 

P.P'.q 

where A is the macroscopic quantization area and C/ q is the Fourier transform of the exciton-exciton pair repulsion 
potential. 

The Hamiltonian of photons in a semiconductor microcavity is given by [24| : 

Hph = 22 e p h[P)a^a P , (5) 
p 

where a P and Sp are the photonic creation and annihilation Bose operators, and e p h{P) is the cavity photon energy 
dispersion. 

Follow Ref. [25] the Hamiltonian of the harmonic exciton-photon coupling can be written as 



3 



H ex - ph = hrt R ^2 ap^p + h.c, (6) 
p 

where the exciton-photon coupling energy represented by the Rabi splitting constant Kin is obtained below. 

Let us present the detailed consideration and analysis of each term in (JXJ) : firstly, the excitonc Hamiltonian H ex that 
describes the formation excitons in gapped graphene, secondly, the Hamiltonian H p h that describes the microcavity 
photons and lastly, the Hamiltonian H ex — P h responsible for the exciton-photon coupling within the microcavity. 



III. EXCITONIC HAMILTONIAN 

In this Section we present the excitonic Hamiltonian H ex in Eq. (flj that consist from two terms: the Hamiltonian 
that describes the formation of the gas of non-interacting excitons in gapped graphene and Hamiltonian responsible 
for the exciton-exciton interaction that we assume is strong enough to be neglected. 



A. An exciton in a graphene layer in the presence of the gap 

The Hamiltonian of 2D excitons in the gapped graphene is given by Eq. @. As the first step we analyze the 
Hamiltonian of non-interacting excitons . As it follows from Eq. @ this Hamiltonian is determined by the 
energy- momentum despersion of non- interacting excitons e ex (P). Therefore, we need to find out the energy- momentum 
despersion e ex (P) of the electron and hole that are bound via an electromagnetic attraction by solving the two-body 
problem in a gapped graphene layer. 

We consider an electron and a hole located in a single graphene sheet and introduce the gap parameter S. We 
assume that our exciton is formed by the electron and hole located in one graphene sheet. The gap parameter i5 is 
the consequence of adatoms on the graphene sheets (e.g, by hydrogen, oxygen or other non-carbon atoms Q) which 
create a one-particle potential. 

We use the coordinate vectors of the electron and the hole ri and r 2 . Each honeycomb lattice is characterized by 
the coordinates (ij, 1) on sublattice A and (rj,2) on sub-lattice B. Then the two-particle wave function, describing 
two particles in the same graphene sheet, reads ^(rx, S\; t%, 82)- This wave function can also be understood as a 
four-component spinor, where the spinor components refer to the four possible values of the sublattice indices si, s 2 : 



*(ri,si;r 2 ,s 2 ) 



( 0Qa(ri,r 2 ) 
(f>ab{ri,r 2 ) 

V 0bft(ri,r 2 ) 



(7) 



The system of two interacting particles located in the same graphene sheet by using the coordinates r, which is 
the projection of the difference vector between between an electron and a hole ri — r 2 on the plane parallel to the 
graphene sheet. Then the Hamiltonian of the two-body problem with a broken sublattice symmetry can be described 
by the Hamiltonian 
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V(r)J 



(8) 



where V(r) is the electron-hole electromagnetic interaction as function of r, which is the projection of the vector 
rx — r 2 on the plane parallel to the graphene sheet. In Eq. flSJ), dj — hvp(—id Xj — d Vj ), d\ — Hvjp(—id x , + d. 



J Vj> 



where d Xj — d/dxj, d Vj = d/dyj, j = 1,2, and x and y are components of the vectors ri and r 2 that represent 
the coordinates of the electron and hole, correspondingly, and Vp — y/3at/(2h) is the Fermi velocity of electrons in 
graphene, where a = 2.566 A is a lattice constant and t ps 2.71 eV is the overlap integral between the nearest carbon 
atoms [26j |. In Eq. ([8]) we take into account the renormalization of the electron- hole distance in the electron- hole 
Coulomb attraction due to the non-locality of electron and hole wave function assuming the following model for 
the electron-hole attraction V(r) — — e 2 /(47reoe(r 2 +^o) 1 ^ 2 )j where ro is the renormalization parameter which will be 
estimated, eo = 8.85 x 1CP 12 C 2 /(Nm 2 ), e is the electron charge, e is the dielectric constant of the material surrounding 
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graphene sheet. It should be mentioned that the main contribution to the polariton mass is the cavity photon mass 
rather than exciton mass, and, therefore, we use our model to estimate the exciton mass roughly by the order of 
magnitude. Assuming r <C Tq, we expand V(r) in Taylor series as V(r) — —Vq + jr 2 , where Vq = e 2 / (Aireearo) and 
7 = e 2 /(87T£ £^). 

Now we have to solve the eigenvalue problem of the Hamiltonian H^B = e^> for the energy of the exciton e. The 
eigenfunction depends on the coordinates of both particles, namely (ri,ra). To separate the relative motion of the 
electron and hole we use the following ansatz for the wave function 

^(R,r) = e ip - R /^-(r) , (9) 

where P is momentum, and following to Ref. [2l] for a generalized center of mass coordinate R and relative coordinate 
r we have 



R = ari + f3r 2 , r = n r 2 , (10) 

with the parameters 

e ~ 26 /3 = ^- (ID 



2e ' 2e 

The procedure given in Ref. l2ll can be applied to the eigenvalue problem (H + Vo)^ = (e+ Vb)* when we assume that 
both relative and center-of-mass kinetic energies, as well as the harmonic term in the potential energy V(r) + Vo = jr 2 
are small in comparison to the gap energy 26. Starting from the wave function in Eq. ([7]) and by using the same 
procedure as we applied in Ref. l2l| we obtain for the spinor component 4> aa the equation 

(v F P) 2 , 2e(hv F ) 2 S7 2 \ 

V ( r ) ~T? THT § aa = 6< Paa ■ (12) 



2e w (e 2 - A5 2 ) 

Let us rewrite Eq. (1121) in more convenient form 

(-J-V 2 . + 7 r 2 ) <P aa = T (f>aa , (13) 

where T is given by 



^_ 2(e + V )(hv F f 

((e + V ) 2 -4S 2 ) [ ' 



and Tq is given by 



T °='- +v °-wfw)' < 15 ' 

Eq. (|13p describes a two-dimensional isotropic harmonic oscillator, whose solutions are given by the condition 

yr = 2iV^ (16) 

with N = 2ni + n 2 + 1 and quantum numbers n\ = 0,1,2,3,..., n 2 = 0, ±1, ±2, ±3, . . . , ±ni. We will focus 
subsequently on the analysis of the ground state corresponding to N = 1. 

We solve Eq. (|16|) for e, then expand e up to the second order in P (i.e. for vpP <C 6) and obtain for the exciton 
dispersion in Eq. (|3]) 



e«(P) = -Vo + 2fyi+£ + |i, (17) 
where C = -f(hvp) 2 and the effective exciton mass obtained under assumption that C « e(e 2 - 45 2 ) is given by 



25_ 4 /„ C 

4< 
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We note the exciton effective mass M. increase when the gap 25 increases. The exciton energy spectrum increases 
with increasing gap 25. Since annihilation of the exciton results in the emission of a photon whose energy is that of 
the gap, the renormalization parameter tq must satisfy the condition e ex (0) = 25. Therefore, 

e 2 1 / 

r ° = -^ x+7x\ La 2 2 2 + 2h2v F ■ (19) 

32ne e5 45 y 647r 2 £ 2 ,£ 2 * 

Then the exciton radius can be obtained from the wavefunction of the 2D harmonic oscillator and reads p — 
l/2y / J r (e(P = 0))/7. It is easy to show that the exciton radius p = 1.5r , where we assumed that the graphene 
layer was surrounded by GaAs with the dielectric constant e = 13. 

We note the exciton effective mass Ai increases when the gap parameter 5 increases. The exciton energy spectrum 
at the same quantum number N increases with the increase of the gap 5 at small momenta P. The exciton energy 
spectrum at the same gap 5 decreases with the increase of the quantum number AT. 



B. Exciton-exciton interaction 



Here we analyze the Hamiltonian of exciton-exciton interaction H ex _ ex in graphene in the presence of gap given 
by Eq. (|4]), which contributes to the exciton Hamiltonian H ex given by Eq. ([2]). As discussed in Refs. [2?] and for a 
dilute exciton gas, the excitons can be treated as bosons with a repulsive contact interaction. For small wave vectors 
q <C p^ 1 the pairwise exciton-exciton repulsion can be approximated as a contact potential J7 q ~ U = 3e 2 pj \2i:eoe). 
This approximation for the exciton-exciton repulsion is applicable, because resonantly excited excitons have very small 
wave vectors [25j . Another reason for the validity of this approximation is that the exciton gas is assumed to be very 
dilute and the average distance between excitons r s ~ (ra) -1 / 2 3> p, which implies the characteristic wavenumber 
q ~ rj 1 <C p~ x ■ A much smaller contribution to the exciton-exciton interaction is also given by band-filling saturation 
effects 29], which are neglected here. 

Thus, since U is directly proportional to p, and p is inversely proportional to 5, U is inversely proportional to 5. 
Therefore, we can conclude that the exciton-exciton interaction U decreases when 5 increases. 



IV. MICROCAVITY PHOTONS 



The Hamiltonian of photons in a semiconductor microcavity is determined by the cavity photon energy dispersion 



e p h(P). According to Ref. \24 this dispersion is defined as e P h(P) — (c/n)y P 2 + h 2 K 2 L c 2 , where n = *Je is the cavity 
effective refractive index that is given by the dielectric constant of the cavity, c is the speed of light in vacuum and 
Lq is the length of the microcavity. 

Embedding graphene in an optical microcavity can lead to the formation polaritons, when the excitons couple to 
the cavity photons. In such a case the microcavity consists of two mirrors parallel to each other and a graphene sheet 
placed in between. Then the photons are confined in the direction perpendicular to the mirrors, but move freely in 
the two directions parallel to the mirrors. The length of the microcavity is chosen as 

Lc = 2M (20) 

with the resonance condition that the photonic and excitonic branches agree at P = 0, i.e. for e ex (0) — e p /i(0). 
This resonance condition can be achieved either by controlling the dispersion of excitons e ex (P) or by choosing the 
appropriate length Lc of the microcavity. 



V. EXCITON-PHOTON INTERACTION 



We derive the Hamiltonian of the harmonic exciton-photon coupling ([6|), which contributes to the total Hamiltonian 
of the system H to t given by Eq. ([T]). In this Hamiltonian exciton-photon coupling energy is represented by the Rabi 
splitting constant KIr. Neglecting anharmonic terms for the exciton-photon coupling, the Rabi splitting constant 
VLr can be estimated quasiclassically as 



I Ml, 



/ H int 



(21) 
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where Hint is the Hamiltonian of the electron-photon interaction, |/) and \i) are the final and initial states of the 
system, correspondingly. The initial state \i) corresponds to the filled by electrons valence band and empty conduction 
band, while the final state |/) corresponds to one hole in the valence band and one electron in the conduction band. 
The eigenfunctions and eigenenergies of an electron in graphene in the presence of gap are given in Appendix |XJ For 
graphene this interaction is determined by Dirac electron Hamiltonian as 



tlint — CT ' Ap/,0 — — a ■ hiphQ 

C IUJ 



(22) 



where a — {a x ,a v ) are Pauli matrices, A p ^o is the electromagnetic vector potential of a single cavity photon, and 

1 /2 

Epho = (2nfrw / (eW)) is the magnitude of electric field corresponding to a single cavity photon with frequency uj 
with the microcavity volume W . 

In Eq. ([2~Tj) the initial \i) and final |/) electron states are defined as 

\i) = n^io^io^, 



l/> = &t|i) 



(23) 



In Eq. (|23|) . Cq ^ is the Fermi creation operator of the electron in the valence band with the wavevector q, |0) c 
denotes the wavefunction of the vacuum in the conduction band, riqCq'^IO)-" corresponds to the completely filled 
valence band, 6 q is the exciton creation operator with the electron in the conduction band c and the hole in the 
valence band v. Following Ref. HO, & q and 6 q for this case are defined as 



E 



A(«)t do) 

L q-q'V 



dc)U(v) 

q-q' 



(24) 



where c q and c q are the Fermi creation and annihilation operators of the electron in the conduction band with the 
the wavevector q. 

3j, and (IA4[) into (|23[) and using the electron-photon interaction Hi nt (p2|) . we finally 



Substituting Eqs. (|2H 
obtain from Eq. ([21]) : 



\m f 



evp 



ev F I5 2 -E 2 
2iu 



dx J dy [i/j* cEc (x,y)<7 ■ E pho ip VtE „(x,y)] 

( hv F (q x + iq y ) hv F (q x - iq y ) 



E 2 



Ephox 



5-E 



5 + E 



-iE, 



phOy 



-hv F {q x + iq y ) hv F (q x - iq y ) 



5-E 



5 + E 



After simplification we obtain from Eq. (|25l) 
1 ' ujEVS 2 - E 2 



(qxE p ho x + q y Epho y ) 2 5 2 + (q y E p ho x — qxEphQ y ) 2 E 2 



(25) 



(26) 



Assuming for simplicity electric field corresponding to a single cavity photon mode is directed along the x axis, we 
obtain 



\Ml R \ = 



ev F 



^-J^s 2 + q 2 y E 2 



yj5 2 + h 2 v 2 F q 2 V Weuj 
In the limit q 0, when w = 2S/H, we obtain from Eq. (|27p 



(27) 



hv F e 



WeS 



(28) 
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VI. SUPERFLUIDITY OF MICROCAVITY POLARITONS 



We can diagonalize the linear part of the total Hamiltonian H to t (without the second term on the right-hand side 
of Eq. P])) by applying unitary transformations 25] and obtain (see Appendix [Bj) : 

H = J2 *Lp(P)%,h + Zup{P)ul,u P , (29) 
p p 

where f P , fp and up and u P are the Bose creation and annihilation operators for the lower and upper polaritons, 
respectively. The energy spectra of the low/upper polaritons are given by Eq. (|B3[) . 

Substituting the polaritonic representation of the excitonic and photonic operators (|B1[) into the total Hamilto- 
nian ([1]) , the Hamiltonian of lower polaritons is obtained [25[ : 

H tot =^Lp(P)l f P h + ^ I £ W P , P ,J P+ / P ,_J P l P ,, (30) 

P PP'.q 

where the energy dispersion of the low polaritons £lp(P) is given by Eq. (|B3|) . and the effective polariton-polariton 
interaction W is given by 

Wp.P'.q = ^^P+q^P'^P'-q^P ) (31) 

where U = 3e 2 p/ {2-keqe) according to Sec. IIIIB1 

At small momenta a = 1/2(M~ 1 + (c/h)Lc /fiTr)P 2 /\MIr\ -C 1, the single-particle lower polariton spectrum 
obtained from Eq. (|B3[) . in linear order with respect to the small parameters a, is 

e (P) « UmL^ - \m a \ + ^- , (32) 



where M p is the effective mass of polariton given by 



M P = 2(M- 1 + ^) . (33) 

The mass of polariton M p as a function of S is presented in Fig. [T] According to Fig. [TJ the mass of polariton M p 
increases when S increases, being directly proportional to i5. 

If we take into account only the lower polaritons corresponding to the lower energy and measure energy relative to 
the P = lower polariton energy {c/n)fvKL^j L — \KIr\, the resulting effective Hamiltonian for polaritons has the form 

p2 MO) 

A * = E + ^ E ll+Jl-J^' . ( 34 ) 

P p P,P',q 

where U^g = \U = 3e 2 p/(87re e), since at small momenta \X P \ 2 ps \Cp\ 2 w 1/2. 

In the dilute limit (np 2 <C 1, where n is the 2D polariton density), at zero temperature Bose-Einstein condensation 
(BEC) of polaritons appears in the system, since Hamiltonian of microcavity polaritons Q34[) corresponds to the 
weakly-interacting Bose gas. The Bogoliubov approximation for the dilute weakly-interacting Bose gas of polaritons 
results in the sound spectrum of collective excitations at low momenta [3l|, [32[: e(P) — c$P with the sound velocity 

1/2 

c s = (U^n/Mpj = (3e 2 pn/(8Tre eM p )) 1/2 . 

The dilute polaritons constructed by excitons in graphene in the presence of the gap and microcavity photons when 
gapped graphene is embedded in an optical microcavity form 2D weakly interacting gas of bosons with the pair short- 
range repulsion. So the superfluid- normal phase transition in this system is the Kosterlitz-Thouless transition [33| . 
and the temperature of this transition T c in a two-dimensional microcavity polariton system is determined by the 
equation: 

_ nh 2 n s {T c ) 

c ~ 2k B M p ' [6b) 

where n s (T) is the superfluid density of the polariton system in a microcavity as a function of temperature T, and 
ks is Boltzmann constant. We obtain the superfluid density as n s — n — n n by determining the density of the normal 
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FIG. 1: The dependence of the mass of polariton M p on S. The polariton mass M p is represented in the units of electron mass 




FIG. 2: The dependence of the Kosterlitz-Thouless transition temperature T c on S and the polariton density n. 



component n n when we follow the procedure |31| as a linear response of the total momentum with respect to the 
external velocity: 



n s = n- 



3C(3) sk%T 3 
2ttH 2 c%M p ' 



(36) 



where s = 4 is the spin degeneracy factor. 

Substituting Eq. (|36|) for the density n s of the superfluid component into Eq. (|35j) . we obtain an equation for the 
Kosterlitz-Thouless transition temperature T c . The solution of this equation is 



32 {M p k B T° 



1/3 



27 V Trh 2 r 



32 (M p k B T° 



1/3- 



27 V 7t/i 2 J 



1 - 1 



c 

2V3 



(37) 



where T° is the temperature at which the superfluid density vanishes in the mean- field approximation (i.e., n s (T ( ?) = 
0), 



1 firti 2 nc'tM r} 



6sC(3) 



1/3 



(38) 



The Kosterlitz-Thouless transition temperature T c as a function of 5 and the polariton density n is presented in 
Fig. [2] It can be seen that T c increases, when S decreases and the polariton density n increases. Fig. [3] presents the 
Kosterlitz-Thouless transition temperature T c as a function of 5 at the different fixed polariton densities n. According 
to Fig. [31 at the same polariton density n, T c decreases when 5 increases, and at the same 5, and T c is higher for 
higher n. 
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FIG. 3: The dependence of the Kosterlitz-Thouless transition temperature T c on S. 

VII. DISCUSSION 

According to Eq. (|33|) . the effective polariton mass M p is mostly determined by the size of the microcavity, which 
depends on the gap 8 according to Eq. (f2"Uf when the exciton and microcavity photon branches of the spectrum are 
in resonance at zero momentum. The effective polariton mass M p is proportional to the gap S. The dependence 
of the superfluid density n s on the gap 6 is determined also by dependence of the sound velocity cs on the exciton 
radius p {eg ~ p 1 ^ 2 )- Since p ~ <5 _1 , we have cs ~ <5 -1 / 2 , and the normal component density n n ~ S 2 according to 
Eq. (I36|) . The gap dependence of the sound velocity is caused by the gap dependence of exciton-exciton interaction 
described in Sec. 1111 Bl Therefore, we conclude that the superfluid density n s and Kosterlitz-Thouless temperature T c 
are decreasing functions of the gap 6 in graphene and increasing functions of the exciton density n. Besides, it follows 
from Eq. (|28l) . that Rabi splitting is decreasing function of the gap as MIr ~ <5 -1 / 2 . 

The advantage of observing the superfluidity and BEC of polaritons formed by gapped graphene excitons and 
microcavity photons in comparison with these formed by quantum well excitons and microcavity photons is based on 
the fact that the superfluidity and BEC of polaritons formed by gapped graphene excitons can be controlled by the 
gap which depends on doping. The Rabi splitting related to the creation of an exciton in a graphene layer also can 
controlled by the doping dependent gap. 

In conclusion, we propose the superfluidity of 2D exciton polaritons formed by the gapped graphene excitons and 
microcavity photons, when the gapped graphene layer is embedded in an optical microcavity. The effective exciton 
polariton mass is calculated as a function of the gap energy in the graphene layers. We demonstrate the effective 
exciton polariton mass increases when the gap increases and it is directly proportional to S. We show that the 
superfluid density n s and the Kosterlitz-Thouless temperature T c increases with the rise of the excitonic density n 
and decreases with the rise of the gap due to 5 dependence of the sound velocity of collective excitations (cs ~ 
and therefore, could be controlled by n and S. We demonstrate that the Rabi splitting related to the creation of an 
exciton in a graphene layer also depends on the gap energy. 
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Appendix A: The eigenfunctions and eigenenergies of an electron in graphene in the presence of gap 

We consider two fermions and ignore their interaction (i.e. V(r) = 0). The electrons in the conduction band, 
described by the spinor wave function tp C E(x,y), and the holes in the valence band, described by the spinor wave 
function ip V E'(x,y), arc solutions of the eigenvalue equations 



H_s?PcE = E%j) cE , HsTpvE' = —E'ip vE > 



(Al) 
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of the Dirac-Weyl Hamiltonian 



S Hvp(d x — id y ) 

H s = I hv F (d x +id y ) -6 1 . (A2) 



In the presence of a gap 25, these solutions are 



, jS-E exp (i (q x x + g „y)) / 

M*,v) = \I^e { 1 J (A3) 



*.(.,v> - ^2^^ig±M)) ( ^fe^ , , ( A4) 



where i? = -\/<5 2 + h 2 v 2 F q 2 is the energy of the electron or the hole, and L x and L y are lengths in a; and y direction, 
correspondingly. This allows us to construct the four components of the spinor in Eq. (J7J from the solutions of Eq. 
Iffrljl as 

0jk(ri,r 2 ) = ip C E,j(ri)tpvE',k( r 2) ti = a,b; k = a,b) (A5) 

which solves the eigenvalue equation for two non-interacting particles with the Hamiltonian H in Eq. @ when 
V(r) = 0: "Ho^ = (E — E')*ff, where ^0 is the Hamiltonian of two non-interacting particles. 

Appendix B: The representation of polariton operators for the Hamiltonian of the exciton-photon system in 

a microcavity 

We express the exciton and microcavity photon operators in terms of polariton operators. The exciton and photon 
operators are defined as [25| 

bp = Xplp — Cpup , dp = Cplp + Xpup , (Bl) 
where lp and up are lower and upper polariton Bose operators, respectively, Xp and Cp are 



K e LP (P) - e ph (P)J J V V m R 

and the energy dispersion of the low/upper polaritons are 

_ e ph {P) + e ex (P) 
£lp/up(^) - ^ 



T ^y/ (e ph (P) - e ex (P))* + 4\m R \* . (B3) 
We note that |^p| 2 and \Cp\ 2 = 1 — |^p| 2 represent the exciton and cavity photon fractions in the lower polariton. 
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